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Recent developments in the uses of QCD spectral sum rules (QSSR) for heavy flavours 

are summarized and updated. QSSR results are compared with the existing data and 
with the ones from alternative approaches. 

1 Introduction 

We have been living with QCD spectral sum rules (QSSR) (or QCD sum rules, or ITEP 
sum rules, or hadronic sum rules...) for 15 years, within the impressive ability of the 
method for describing the complex phenomena of hadronic physics with the few universal 
"fundamental" parameters of the QCD Lagrangian (QCD coupling a s , quark masses and 
vacuum condensates built from the quarks and/or gluon fields, which parametrize the non- 
perturbative phenomena). The approach might be very close to the lattice calculations as 
it also uses the first principles of QCD, but unlike the case of the lattice, which is based on 
sophisticated numerical simulations, QSSR is quite simple as it is a semi-analytic approach 
based on a semiperturbative expansion and Feynman graph techniques implemented in 
an Operator Product Expansion (OPE), where the condensates contribute as higher- 
dimension operators. The QCD information is transmitted to the data via a dispersion 
relation obeyed by the hadronic correlators, in such a way that in this approach, one 
can really control and in some sense localize the origin of the numbers obtained from the 
analysis. With this simplicity, QSSR can describe in an elegant way the complexity of 
the hadron phenomena, without waiting for a complete understanding of the confinement 
problem. 

One can fairly say that QCD spectral sum rules already started, before QCD, at the time 
of current algebra, in 1960, when different ad hoc super convergence sum rules, especially 
the Weinberg and Das-Mathur-Okubo sum rules, were proposed but they came under 
control only with the advent of QCD [|TJ. However, the main flow comes from the classic 
paper of Shifman-Vainshtein-Zakharov 0] (hereafter referred to as SVZ), which goes 
beyond the naive perturbation theory thanks to the inclusion of the vacuum condensate 
effects in the OPE (more details and more complete discussions of QSSR and its various 
applications to hadron physics can be found, for instance, in ||). 

1 This is an updated version of the talks given by the author at the XXIXth Rencontre de Moriond, 
Meribel (1994), CERN-TH.7277/94 (1994) and at the QCD94 Conference, Montpellier (1994), CERN- 
TH. 7444/94 (1994). An extended version of this review will be published in Recent Developements of 
QCD spectral sum rules by World Scientific Company. 
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In this talk, I shall present aspects of QSSR in the analysis of the properties of heavy 
flavours. As I am limited in space-time, I cannot cover in detail here all QSSR applications 
to the heavy-quark physics. I will only focus on the following topics, which I think are 
important in the development of the understanding of the heavy-quark properties in 
connection with the progress done recently in the heavy quark effective theory (HQET) 
(or infinite mass effective theory (IMET)) and in lattice calculations: 

- the heavy-quark-mass values, -the meson-quark mass difference and the heavy quark 
kinetic energy, 

- the pseudoscalar decay constants, - the bag constants and the CP- violation parameters, 

- the heavy to light exclusive decays, 

- slope of the Isgur-Wise (IW) function and determination of V c b, 

- properties of the hybrids and S c -like hadrons. 



2 QCD spectral sum rules 

In order to illustrate the QSSR method in a pedagogical way, let us consider the two-point 
correlator: 

nf = i J d'x (0\TJ£(x) (J 6 ») f |o> 

= -(^V-gV)n 6 (g 2 ,M 6 2 ), (1) 

where J£{x) = b^^b{x) is the local vector current of the 6-quark. The correlator obeys 
the well-known Kallen-Lehmann dispersion relation: 

poo pjf 1 

Ml 2 , Ml) = / _ - lmU b (t) + (2) 

Jam* t — q l — it tt 

where ... represent subtraction points. This sum rule expresses in a clear way the duality 
between the spectral function Im II&(t), which can be measured experimentally, as here 
it is related to the e + e~ into T-like states total cross-section, while Ub(q 2 ,M£) can be 
calculated directly in QCD, even at q 2 = 0, thanks to the fact that Ml - q 2 > A 2 . The 
QSSR is an improvement on the previous dispersion relation. 

On the QCD side, such an improvement is achieved by adding to the usual perturbative 
expression of the correlator, the non-perturbative contributions as parametrized by the 
vacuum condensates of higher and higher dimensions in the OPE || : 

Mi 2 , mS) * Y. 1 c/3 

D=o,2,4,... {M b — q z J 
. £ C^(q 2 ,Mlu)(0(u)), 

dimO=D 

(3) 

where v is an arbitrary scale that separates the long- and short-distance dynamics; 

are the Wilson coefficients calculable in perturbative QCD by means of Feynman diagrams 
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techniques: D = corresponds to the case of the naive perturbative contribution; (O) 
are the non-perturbative condensates built from the quarks or/and gluon fields. For 
D — 4, the condensates that can be formed are the quark M^ip) and gluon (a s G 2 ) 
ones; for D = 5, one can have the mixed quark-gluon condensate (ipa^X 0, /2G^ u ip) , while 
for D = 6 one has, for instance, the triple gluon gf abc (G a G b G c ) and the four-quark 
a s {i()Ti%p%l)T2ip)-, where are generic notations for any Dirac and colour matrices. The 
validity of this expansion has been understood formally, using renormalon techniques 
(IR renormalon ambiguity is absorbed into the definitions of the condensates) Q and 
by building renormalization-invariant combinations of the condensates (Appendix of 
and references therein). The SVZ expansion is phenomenologically confirmed from the 
unexpected accurate determination of the QCD coupling a s and from the measurement 
of the QCD condensates from semi-inclusive tau decays and spectral moments 0, ||. In 
the present case of heavy-heavy correlators the OPE is much simpler, as one can show 
0-0 that the heavy-quark condensate effects can be included into those of the gluon 
condensates, so that, up to D < 6, only the (a s G 2 ) and g{G 3 ) condensates appear in the 
OPE. Indeed, SVZ have, originally, exploited this feature for their first estimate of the 
gluon condensate value, though the validity of their result has been criticized later on 



For the phenomeno logical side, the improvement comes from the uses of either a finite 
number of derivatives and finite values of q 2 (moment sum rules): 



M [n) = 



1 d n Ii b {q 



2\ 



n\ {dq 2 ) n 
dt 1 

AM? t n+l 7T 



Imn 6 (t), (4) 



or an infinite number of derivatives and infinite values of q 2 , but keeping their ratio fixed 
as t = n/q 2 (Laplace or Borel or exponential sum rules): 

poo I 

£(t,M 2 )= dt exp(-tr) - lmU b (t). (5) 

J4A1 2 ' 7T 

b 

There also exist non-relativistic versions of these two sum rules, which are convenient 
quantities to work with in the large-quark- mass limit. In these cases, one introduces 
non-relativistic variables E and t^: 

t = {E + M b ) 2 and t n = 4M b r. (6) 

In the previous sum rules, the gain comes from the weight factors, which enhance the 
contribution of the lowest ground-state meson to the spectral integral. Therefore, the 
simple duality ansatz parametrization: 

"one resonance"5(t - M 2 R ) + "QCD continuum"0(t - t c ), (7) 

of the spectral function, gives a very good description of the spectral integral, where the 
resonance enters via its coupling to the quark current. In the case of the T, this coupling 
can be defined as: 

(Ol&y&IT) = V2^. (8) 
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The previous feature has been tested in the light-quark channel from the e + e~ — > I = 1 
hadron data and in the heavy-quark ones from the e + e~ — > T or ip data, within a good 
accuracy. To the previous sum rules, one can also add the ratios: 

^ (B) ^ and n ^~i^ (») 

and their finite energy sum rule (FESR) variants, in order to fix the squared mass of the 
ground state. In principle, the pairs (n,t c ), (r,t c ) are free external parameters in the 
analysis, so that the optimal result should be insensitive to their variations. Stability 
criteria, which are equivalent to the variational method, state that the best results should 
be obtained at the minimas or at the inflexion points in n or r, while stability in t c is 
useful to control the sensitivity of the result in the changes of t c values. To these stability 
criteria are added constraints from local duality FESRs, which correlate the t c value to 
those of the ground state mass and coupling |12| . Stability criteria have also been tested 



in models such as the harmonic oscillator, where the exact and approximate solutions 



are known fLOfl . The most conservative optimization criteria, which include various types 
of optimizations in the literature, are the following: the optimal result is obtained in 
the region, starting at the beginning of r/n stability (this corresponds in most of the 
cases to the so-called plateau often discussed in the literature, but in my opinion, the 
interpretation of this nice plateau as a sign of a good continuum model is not sufficient, 
in the sense that the flatness of the curve extends in the uninteresting high-energy region 
where the properties of the ground state are lost), until the beginning of the t c stability, 
where the value of t c more or less corresponds in some cases to the one fixed by FESR 
duality constraints. The earlier sum rule window introduced by SVZ, stating that the 
optimal result should be in the region where both the non-perturbative and continuum 
contributions are small, is included in the previous region. Indeed, at the stability point, 
we have an equilibrium between the continuum and non-perturbative contributions, which 
are both small, while the OPE is still convergent at this point. 



3 The heavy- quark- mass values 



3.1 The running masses 



Here, we will summarize the recent results obtained in [ Hf , where an improvement and 
an update of the existing results have been done, with the emphasis that the apparent 
discrepancy encountered in the literature is mainly due to the different values of a s used 



by various authors. Using the world average value a s (M z ) = 0.118±0.006 |15|, [L6|, [PTHand 
a conservative value (a^G 2 ) = (0.06 ± 0.03) GeV 4 M, H, the first direct determination of 



the running mass to two loops in the MS'-scheme, from the \l/ and T systems, is [H] : 



m c (Mf T2 ) = (1.23±°;° 2 ± 0.03) GeV 

m b (M b PT2 ) = (4.23l°;|* ± 0.02) GeV, (10) 
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where the errors are respectively due to a s and to the gluon condensate. Using the 
previous result in (10) and the expression of the running mass to two-loops [I]. |3|]: 




m Q {v)=m Q -ft^ x i + £ £_£ ^ , (11) 



in terms of the invariant mass ttiq, one can extract the running mass at another scale; 
71 = 2 and 72 = 101/12 - 5n//18, ft = -11/2 + n//3, ft = -51/4 = 19n//12 are the 
mass anomalous dimensions and the ftfunction in the MS'-scheme. Then, one obtains at 
1 GeV: 

m c (l GeV) = (1.46lg;[j§ ± 0-03) GeV 

fn b (l GeV) = (6.37±g:|| ± 0.07) GeV, (12) 

By combining the previous value of the running 6-quark mass with the s-quark one eval- 
uated at 1 GeV, which we take in the range: m s {l GeV)= 150-230 MeV [18|, [19[, one 
obtains the scale-independent ratio: 



m b /m s ~ 33.5 ± 7.6, (13) 
a result of great interest for model-building and SUSYGUT-phenomenology. 



3.2 The pole masses 

One can transform the results on the running masses into the perturbative pole masses 
by using the perturbative relation [18]: 



M Q (v) = m Q (u 



(14) 



where the constant term of the (a s /7r) 2 is known to be: K b ~ 12.4, K c ~ 13.3 |20|.Then, 
we obtain, to two-loop accuracy: 



M, 



PT2 



PT2 



(1.42 ± 0.03) GeV 
(4.62 ± 0.02) GeV. 



(15) 



It is informative to compare these values with the ones from the pole masses from non- 
relativistic sum rules to two loops: 



M NR 
rNR 



:i.45l°"^ ± 0.03) GeV 



-0.03 
-0.02 



M£ n = (4.69lg;g? ± 0.02) GeV. 
A similar comparison can be done at three-loop accuracy. One obtains: 



(16) 



M, 



PT3 



PT3 



(1.62 ±0.07 ±0.03) GeV 
(4.87 ±0.05 ±0.02) GeV, 



(17) 



to be compared with the dressed mass: 



(4.94 ±0.10 ±0.03) GeV, 
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obtained from a non-relativistic Balmer formula based on a bb Coulomb potential and 
including higher order a s -corrections [ETfj . One can remark that the radiative a 2 s correction 



is large and causes a positive shift of about 250 MeV on the value of the pole mass M b . One 
can also remark that at the two and three loop-accuracies, the mass- difference between the 
relativistic and non-relativistic pole masses is about 70 MeV. The interpretation of this 
mass-difference is not quite well understood. If one has in mind that the non-relativistic 
pole mass contains a non-perturbative piece due to Coulombic interactions, which can be 
of the same origin as the one induced by the truncation of the perturbative series at large 
order, then one can consider this value as a phenomenological estimate of the renormalon 
contribution, which is comparable in strength with the estimate of about 100-133 MeV 
from the summation of higher order corrections of large order perturbation theory |22| . 
An extension of the previous analysis of the and T-systems to the case of the B and 
B* mesons leads to the value: 

M b PT2 = (4.63 ± 0.08) GeV, (19) 
in good agreement with the previous results, but less accurate. 



3.3 The b and c pole-mass-difference 

One can also use the previous results, in order to deduce the mass-difference between the 
b and c (non)-relativistic pole masses: 

M b (M b ) - M C (M C ) = (3.22 ± 0.03) GeV, (20) 

in good agreement (within the errors) with potential model expectations |23], |16| . A direct 



comparison of this mass-difference with the one from the analysis of the inclusive -B-decays 
needs however a better understanding of the mass definition and of the value of the scale 
entering into these decay-processes. Indeed, if one chooses to evaluate these pole masses 
at the scale v = M b , which can be a natural scale for this process, one obtains to two-loop 
accuracy: 

M e (v = M b ) = (1.08 ± 0.04) GeV, (21) 
which leads to the mass-difference: 

M b - M c \ v=Mb = (3.54 ± 0.05) GeV, (22) 
in good agreement with the one extracted from the analysis of the inclusive -B-decays 



4 The meson- quark- mass gap and the heavy-quark- 
kinetic energy 

The meson-quark mass gap A is in important input in HQET (IMET) approach. It can 
be defined as [[25[] []: 

M B = M 6 + A--|-(^ + 3£), (23) 

^M b 

2 We are aware of the fact that in the lattice calculations, A defined in this way can be affected by 
renormalons 
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where: 



K 



2M 



B 



-{B{v)\K\B{v)) and 



6M, 



■{B(v)\S\B(v)) 



(24) 



correspond respectively to the matrix elements of the kinetic and of the chromomagnetic 
operators: 



/C = h{iDfh and 



S = -ha^F^h, 



(25) 



where h is the heavy quark field and F^ v the electric field tensor. The estimate of A from 
HQET-sum rules leads to 



A ~ (0.52 - 0.70) GeV, 



(26) 



in good agreement with the previous results 0, p£i|1 , though less accurate as we have 
taken a larger range of variation for the continuum energy. An anologous sum rule in the 



full QCD theory leads to [ 30 



A ~ (0.6 - 0.80) GeV, 



which combined together leads to the intersecting range of values 

A ~ (0.65 ±0.05) GeV. 

The sum rule estimate of the kinetic energy gives P7fl : 

K ~ -(0.5 ±0.2) GeV 2 



(27) 



(28) 



(29) 



where the large error, compared with the previous result of pi] , is due to the absence of 
the stability point with respect to the variation of the continuum energy threshold. By 
combining the previous estimates with the one of the chromomagnetic energy: 



S^ i (Ml,-M 2 ), 



(30) 



one deduces the value of the pole mass to two-loop accuracy: 

M b = (4.61 ± 0.05) GeV, (31) 
in good agreement with the previous values from the sum rules in the full theory and 



(within the errors) with the earlier HQET results of [28 



5 The pseudoscalar decay constants 

5.1 Estimate of the decay constants 

The decay constants fp of a pseudoscalar meson P are defined as: 

K + M Q )(0\q(i lB )Q\P) = V2M 2 P f P , (32) 
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where in this normalization f n = 93.3 MeV. A rigorous upper bound on these couplings 
can be derived from the second-lowest super convergent moment: 



(33) 



2! (dq 2 ) 2 

where ^ 5 is the two-point correlator associated to the pseudoscalar current. Using the 



positivity of the higher-state contributions to the spectral function, one can deduce [32 



M P [ m a | 

^17 { 1 + 3 m5 + ' 751<s ' + -}' (34) 

where one should not misinterpret the mass- dependence in this expression compared to 
the one expected from heavy-quark symmetry. Applying this result to the D meson, one 
obtains: 

Sd < 2.14/;. (35) 

Although presumably quite weak, this bound, when combined with the recent determi- 
nation to two loops |33| : 

(1.15 ± 0.04) U, (36) 

JD 

implies 

/ Ds < (2.46 ±0.09)/,, (37) 

which is useful for a comparison with the recent measurement of fo s from WA75: fo s — 
(1.76 ± 0.52)/, and from CLEO: f Ds ~ (2.61 ± 0.49)/,. One cannot push, however, 
the uses of the moments to higher n values in this D channel, in order to minimize the 
continuum contribution to the sum rule with the aim to derive an estimate of the decay 
constant because the QCD series will not converge at higher n values. In the D channel, 
the most appropriate sum rule is the Laplace sum rule. The results from different groups 
are consistent for a given value of the c-quark mass. Using the table in |33| and the value 



of the perturbative pole mass obtained previously, one obtains to two loops: 

f D ~ (1.35 ± 0.04 ±0.06)/, =>- f Ds ~ (1.55 ±0.10)/,. (38) 

For the B meson, one can either work with the Laplace, the moments or their non- 
relativistic variants. Given the previous value of M^, these different methods give con- 
sistent values of Jb, though the one from the non-relativistic sum rule is very inaccurate 
due to the huge effect of the radiative corrections in this method. The best value comes 



from the Laplace sum rule; from the table in ||33|| , one obtains: 

f B ~ (1.49 ±0.06 ±0.05)/,, (39) 

while |33|1 : 



^-1.16 ±0.04, (40) 

JB 

where the most accurate estimate comes from the "relativistic" Laplace sum rule. The 
apparent disagreement among different existing QSSR numerical results in the literature 
is not essentially due to the choice of the continuum threshold as misleadingly claimed in 
the literature but is mainly due to the different values of the quark masses used because 
the decay constants are very sensitive to that quantity as shown explicitly in 
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5.2 Static limit and 1/M-corrections to f B 



One could notice, since the first result fs — jo of |34| , a large violation of the scaling 
law expected from heavy-quark symmetry. This is due to the large 1 / M;,-correction found 
from the HQET sum rule |29| and from the one in full QCD |30, [35]. Using the estimate 



of the decay constant in the static limit [27 



/^~(1.98±0.31)A 



(41) 



and the previous estimates of fs and fr> in the full theory, the quark-mass dependence of 
the decay constant can be parametrized as: 



f B jM b ~ (0.33 ± 0.06) GeV 3/2 a s 1//3 \ \ 1 



2 a, 



A B 
+ 



3 7i M b ' Ml j ' 



by including the quadratic mass corrections, where: 

A « 1.1 GeV and B « 0.7 GeV 2 , 
while a linear parametrization leads to: 

A~ (0.6 ±0.1) GeV, 



(42) 



(43) 



(44) 



in accordance with previous findings |2^, |35] and with the lattice results |36| . One can 
qualitatively compare this result with the one obtained from the analytic expression of the 
moment or from the semilocal duality sum rule, which leads to the interpolating formula 

m 

Ell 2 1/fl . /M*\ 3/2 L 2 a, 3 El n 2 (uu) 



+ 



88M fe 2 



2 El 



+ 



(45) 



and gives: 



-(M B -M 6 )~lGeV, 
3 ^ 2 9 



A 

2 ' 

B « - "-{M B - M b y ~ 0.5 GeV* 

in agreement with the previous numerical estimate. 



(46) 



6 The bag constants and the CP-violation parame- 
ters 



6.1 Estimate of the bag constant B B 

The B°-B° mixing is gouverned by the ^-parameter as: 



(47) 



where one can introduce the invariant bag parameter B B as: 



B B = B B {p){a a {v))-^\ (48) 

We have tested the validity of the vacuum saturation B B = 1 of the bag constant, using 
a sum rule analysis of the four-quark two-point correlator to two loops |38[ following 



the leading order work of p9 [. We found that the radiative corrections due to the non- 
factorizable contributions are quite small. Under some physically reasonable assumptions 
for the spectral function, we found that the vacuum saturation estimate is only violated 
by about 15%, giving: 

fi B ~l±0.15. (49) 
By combining this result with the one for f B , we deduce: 

197 ± 18) Mel/, (50) 

if we use the normalization f n = 132 MeV, which is v2 times the one defined in (30), in 
excellent agreement with the present lattice calculations EI 



6.2 Estimate of the bag constant B 



K 



We have also estimated the .B^-parameter associated to the K°-K° mixing, using the four- 
quark two-point correlator as in |40|. Using the Laplace sum rule (LSR) and adopting 



the parametrization of the spectral function in |4(|, we have obtained the conservative 



estimate [O]: 

B K ~ (0.58 ±0.22), (51) 

in good agreement (within the errors) with the FESR result (0.39 ± 0.10) and with ones 
from other approaches |[42|| . However, our central value is slightly higher than the one from 
FESR, where the latter result is mainly due to the effects of the higher radial excitations 
in the FESR analysis which are not under good control. LSR is less sensitive to these 
effects due to the exponential factor which suppresses their relative contributions. One 
can also notice that this result from the two-point function sum rule is more accurate 
than the one from the three-point function which ranges from 0.2 to 1.3, though 

the result of [43] is in good agreement with ours. This inaccuracy can be intuitively 
understood from the relative complexity of the three-point function sum-rule analysis. 



6.3 Estimate of the CP-violation parameters (/?, ij) 

We are now ready to discuss the implications of the previous results for the estimate of 
the CP-violation parameters (p, n) defined in the standard way within the Wolfenstein 



parametrization [[L6|, |45| . Using the previous values of f B , B B and B K , which are all of 
them obtained from a Laplace sum rule analysis, and using the other input used in |57j, 
one obtains the best fit pl|j : 

(p,tj)« (0.09,0.41), (52) 

in very good agreement with the expectation in [Q derived from an alternative method 
(see also |H2[ ^7|). Here, the value of p is very sensitive to the change of B K and f B . 
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7 The heavy to light exclusive decays 



7.1 Introduction and notations 



One can extend the analysis done for the two-point correlator to the more complicated 
case of three-point function, in order to study the form factors related to the B — > K*^ 
and B — > p/rr semileptonic decays. In so doing, one can consider the generic three-point 
function: 

V(p,p',q 2 ) = -Jd A x d 4 y e^-w) (0\TJ L (x)O(0)f B (y)\0), (53) 

where Jl, Jb are the currents of the light and B mesons; O is the weak operator specific 
for each process (penguin for the K*j, weak current for the semileptonic); q = p — p' 3 . 
The vertex obeys the double dispersion relation : 



w> - Q T^hr, Q 5? Iml ^' s '» (54) 



As usual, the QCD part enters in the LHS of the sum rule, while the experimental 
observables can be introduced through the spectral function after the introduction of 
the intermediate states. The improvement of the dispersion relation can be done in 
the way discussed previously for the two-point function. In the case of the heavy to 
light transition, the only possible improvement with a good M& behaviour at large Mj, 
(convergence of the QCD series) is the so-called hybrid sum rule (HSR) corresponding to 
the uses of the moments for the heavy-quark channel and to the Laplace for the light one 

1 r°° ds 



H(n,r') = — / — — / ds' e~ TS lmV(s,s'). (55) 

7T 2 J MS S n+1 JO 



The different form factors entering the previous processes are defined as: 



MP0I«7 M (1 - K)b\B(p)) = (M B + M P )A< - 



A 2 2V 



-e*p'(p + P % + 



Mb + M p M B + M p 

(7v(p')\u^pb\B(p)} = f+(p + p) p + f-(p-p')p, 

(56) 



and: 



(p(p')I^V (^p) ^\B(p)) 



3 It has to be noticed that we shall use here, like in |47|-[p0|, the pseudoscalar current Jp = (m u + 
md)u(ij 5 )d for describing the pion, where the QCD expression of the form factor can be deduced from 
the one in (5l) by taking m c — and by remarking that the additionnal effect due to the light quark 
condensate for B — > ir relative to B — > D vanishes in the sum rule analysis. In the literature |52|, p3| . 
the axial-vector current has been used. However, as it is already well-known in the case of the two-point 
correlator of the axial- vector current, by keeping its q p q v part, (which is similarly done in the case of 
the three-point function) one obtains the contribution from the it plus the A\ mesons but not the 7r 
contribution alone. Though, the A± effect can be numerically small in the sum rule analysis due to its 
higher mass, the mass behaviour of the form factor obtained in this way differs significantly from the one 
where the pseudoscalar current has been used due to the different QCD expressions of the form factor in 
the two cases. 
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* We* WF^ P + 



{e;(M B -M 2 p )-e*q(p + p%}^. 



(57) 



7.2 q 2 and M^-behaviours of the form factors 



We have studied analytically the previous form factors ||48||- ||50|| . We found that they are 



dominated, for M b — > oo, by the effect of the light-quark condensate, which dictates the 
Mb behaviour of the form factors to be typically of the form: 




where Xp is the integral from the perturbative triangle graph, which is constant as 
t' 2 E c / (dd) (t' c and E c are the continuum thresholds of the light and b quarks) for large 
values of Mb. It indicates that at q 2 = and to leading order in 1/M&, all form factors 
behave like \JM\>, although, in most cases, the coefficient of the 1/M 2 term is large. The 
study of the q 2 behaviours of the form factors shows that, with the exception of the A\ 
form factor, their q 2 dependence is only due to the non- leading (in 1/Mj,) perturbative 
graph, so that for M b — > oo, these form factors remain constant from q 2 = to q^ax- The 
resulting M b behaviour at g^ a2 . is the one expected from the heavy quark symmetry. The 
numerical effect of this g 2 -dependence at finite values of Mf, is a polynomial in q 2 (which 
can be resummed), which mimics quite well the usual pole parametrization for a pole 
mass of about 6-7 GeV. The situation for the A\ is drastically different from the other 
ones, as here the Wilson coefficient of the (dd) condensate contains a q 2 dependence with 
a wrong sign and reads [35]: 



Atft-IML-*], (59) 



Sb { M SJ 

which, for q^^ = (Mb — M p ) 2 , gives the expected behaviour: 



M<LJ - t^- (6°) 



It should be noticed that the q 2 dependence of A\ is in complete contradiction with the 
pole behaviour due to its wrong sign. This result explains the numerical analysis of 
]53| . One should notice that a recent phenomenological analysis of the data on the large 



longitudinal polarization observed in B — * K* + \& and a relatively small ratio of the rates 
B — ► K* + over B — > K + ^ can only be simultaneously explained if the A\(q 2 ) 
form factor decreases |j4| as expected from our previous result, while larger choices of 



increasing or/and monotonically form factors fail to explain the data [56|. It is still urgent 
and important to test this anomalous feature of the Ax-form factor from some other data. 
It should be finally noticed that owing to the overall 1/fs factor, all form factors have a 
large 1/M& correction. 
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7.3 Numerical estimate of the form factors and decay rates 



In the numerical analysis, we obtain at q 2 = 0, the value of the B — > K*^ form factor 



48 



F^ p ~ 0.27 ±0.03, 



if" 



1.14 ±0.02, (61) 



which leads to the branching ratio (4.5 ±1.1) x 10~ 5 , in perfect agreement with the CLEO 
data and with the estimate in J57|. One should also notice that, in this case, the coefficient 
of the 1/Mfc correction is very large, which makes dangerous the extrapolation of the c- 
quark results to higher values of the quark mass. This extrapolation is often done in some 
lattice calculations. 

For the semileptonic decays, QSSR give a good determination of the ratios of the form 
factors with the values for the B-decays|[47||: 

A 2 (0) V(0) 



A(o) ^(o) 

Ax(0) 



1.11 ±0.01 



/+(0) " 



~ 1.18 ±0.06 
1.40 ±0.06, (62) 



though their absolute values are inaccurate [47, 53 . This is due to the cancellation of 



systematic errors in the ratios. Combining these results with the "world average" value 
of /+(0) = 0.25 ± 0.02 and the one of F 1 B ^ P (0), one can deduce the rates: 

I\~(4.3±0.7)|K b | 2 xl0 12 s- 1 

r p/ T\ ~ 0.9 ±0.2 (63) 

These results are quite precise and indicate the possibility to reach V u b with a good 
accuracy from the exclusive modes. One should notice here, mainly because of the non- 
pole behaviour of Af , the ratio between the widths into p and into 7r is about 1, while in 
different pole models, it ranges from 3 to 10. Recent data on5-> K(K*) + ^(^') decays 
f)5l favour this result. For the asymmetry, one obtains a large negative value of a, contrary 



to the case of the pole models. 



7.4 SU(3) breaking in B/D — > Klv and determination of V c d/V cs 

We extend the previous analysis for the estimate of the SU(3) breaking in the ratio of 
the form factors: 

R P = /^(0)//^(0), (64) 
where P = B, D. As mentioned before, we use the hybrid moments for the B and the 



double exponential sum rules for the D. The analytic expression of Rp is given in |50 
which leads to the numerical result: 

R B = 1.007 ± 0.020 R D = 1.102 ± 0.007, (65) 
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where one should notice that for M& — > oo, the SU(3) breaking vanishes, while its size 
at finite mass is typically of the same sign and magnitude as the one of fo a or of the 
B — > K*^j discussed before. What is more surprising is the fact that using the previous 
value of Rd with the present value of CLEO data |58|| : 



Br(D- 



7T 



°lv) 



Br(D~ 



K°lu) 



(8.5 ±2.7 ±1.4)%, 



one deduces 



(66) 



(67) 



V cd /V cs = 0.322 ± 0.056, 
Using \Vcd\ = 0.204 ± 0.017 from PDG pf, one then obtains: 

Vcs = 0.63 ± 0.12. (68) 

We can also determine directly the absolute value of the D —> K form factor. We obtain: 

f+^ K (0) ~ 0.80 ±0.16, (69) 



which used into the CLEOII data [16]: 



tf-*(0) 



0.495 ±0.036, 



leads to: 



Vcs = 0.88 ±0.18. 
The average of our two determinations is: 

V„ = 0.71 ±0.10, 



(70) 
(71) 

(72) 



which needs a confirmation of the CLEOII data. One can compare this value with the 
one quoted by PDG94 [jlB]] . We expect that the most reliable result is the lower bound 
derived from Eq. (70) and from fP~* K (0) < 1, which is: 



V„ > 0.62, 



(73) 



while the value V cs ~ 1.01±0.18 quoted there is related to the choice f+^ K (0) ~ 0.70±0.1. 



8 Slope of the Isgur— Wise function and determina- 
tion of Vcb 



Let me now discuss the slope of the Isgur-Wise function. Taron-de Rafael [60] have 
exploited the analyticity of the elastic 6-number form factor F defined as: 

(B(p')\h%\B(b)) = (p + P TF(q 2 ), (74) 



J The old MARKIII data KM would imply a value 0.25 ± 0.15. 
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which is normalized as F(0) = 1 in the large mass limit Mb — Md- Using the positivity 
of the vector spectral function and a mapping in order to get a bound on the slope of F 
outside the physical cut, they obtained a rigorous but weak bound: 

F'(vv' = 1) > -6. (75) 

Including the effects of the T states below BB thresholds by assuming that the TBB 
couplings are of the order of 1, the bound becomes stronger: 

F'(vv' = 1) > -1.5. (76) 

Using QSSR, we can estimate the part of these couplings entering in the elastic form 
factor. We obtain the value of their sum pLf : 



$> T s B ~ 0.34 ±0.02. (77) 

In order to be conservative, we have multiplied the previous estimate by a factor 3 larger. 
We thus obtain the improved bound: 

F'(vv' = 1) > -1.34, (78) 

but the gain over the previous one is not much. Using the relation of the form factor 
with the slope of the Isgur-Wise function, which differs by —16/75 log a s (M b ) one 
can deduce the final bound: 

C'(l) > -1.04. (79) 

However, one can also use the QSSR expression of the Isgur-Wise function from vertex 
sum rules p9| in order to extract the slope analytically. To leading order in 1/M, the 
physical IW function reads: 

C Phys (y = w>) = (_y a |i + ^/ (y ) 

- (dd)r 3 g(y) + {a s G 2 )r A h{y) 
+ g(dGd)r 5 k(y)^, 

(80) 

where r is the Laplace sum rule variable and /, h and k are analytic functions of y. From 
this expression, one can derive the analytic form of the slope ||61|| : 

CphysiV = 1) - -1 + <W + &NP, (81) 

where at the r-stability region: S per t — —Snp — —0.04, which shows the near-cancellation 
of the non-leading corrections. Adding a generous 50% error of 0.02 for the correction 
terms, we finally deduce the leading order result in 1/M: 

C phys (y = 1) - -1 ± 0-02. (82) 

Using this result in different existing model parametrizations, we deduce the value of the 
mixing angle, to leading order in 1/M: 



' cb 



1 48 ns\ x l 2 

■ 1 > x (37.3 ±1.2 ±1.4) x 10" 3 , (83) 



n 
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where the first error comes from the data and the second one from the model-dependence. 



Let us now discuss the effects due to the 1/M corrections. In so doing, we combine the 
predicted value of the form factor^ 0.91 ± 0.03 at y=l, with the one 0.53 ± 0.09 f] from 
the sum rule in the full theory (without a 1/M-expansion) at q 2 = j47|. The model 



dependence of the analysis enters through the concavity of the form factor between these 
two extremal boundaries. We use a linear parametrization: 



C = Co + C'(y-i), 



(84) 



which is also supported by the CLEO data |63| . Then, we can deduce the slope: 

C' ~ -(0.76 ±0.2). (85) 

It indicates that the 1/M correction tends also to decrease the value of which implies 
that, for larger values of y where the data are more accurate, the increase of V c b is weaker 
(+ 3.7%) than the one at y = 1. This leads to the final estimate: 



V, 



cb 



1.48 ps 



n 



1/2 



x (38.8 ± 1.2 ± 1.5 ± 1.5) x 10" 



36) 



where the new last error is induced by the error from the slope, while the model depen- 
dence only brings a relatively small error. Our results for the slope and for V& are in 
good agreement with the new CLEO data f63fl . However, despite its model dependence, 



we expect that the result for V ch is more precise than the one obtained by exploiting the 
value of the Isgur-Wise function at y — 1 [Ml, where the data near this point are quite 



inaccurate. It also shows that the value from the exclusive channels is slightly lower than 



the present result from the inclusive mode p4 |, which is largely affected by the large 
uncertainty in the quark-mass definition and in the heavy quark kinetic energy entering 
into the inclusive process. 



9 Properties of the hybrids and i^-like hadrons 

Let me conclude this talk by shortly discussing the masses of the hybrid QGQ and the 
mass and decays of the 5 c like-hadrons. 



9.1 The hybrids 

Hybrid mesons are interesting because of their exotic quantum numbers. Moreover, it is 
not clear if these states are true resonances or if they only manifest themselves as a wide 
continuum instead. The lowest cGc states appear to be a l + ~ of mass around 4.1 GeV ||. 
The available sum- rule analysis of the 1~ + state is not very conclusive due to the absence 
of stability for this channel. However, the analysis indicates that the spin-one states are in 
the range 4.1-4.7 GeV, in agreement with the predictions from alternative methods |)5 



5 We have taken a compromise value between the ones in 
6 This value is just on top of the CLEO data []63[. 
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Their characteristic decays should occur via the rj' £/(l)-like particle produced together 
with a if) or an rj c . However, the phase-space suppression can be quite important for these 
reactions. The sum rule predicts that the , ++ cGc states are in the range 5-5.7 
GeV, i.e. about 1 GeV above the spin one. Intensive searches of these particles in the 
next r-charm and B factories are an alternative test of our idea about the confinement of 
QCD. 

9.2 The Sc-like hadrons 

We have estimated the S c -meson mass and coupling by combining the results from poten- 
tial models and QSSR ||. We predict the spectra of the 5 c -like hadrons from potential 
models: 

(6255 ± 20) MeV, 
(6330 ± 20) MeV, 
(6.93 ± 0.05) GeV, 
(7.00 ± 0.05) GeV, 
(3.63 ± 0.05) GeV, 

(10.21 ±0.05) GeV, (87) 

which are consistent with, but more precise than, the sum-rule results. The decay constant 
of the B c meson is better determined from QSSR. The average of the sum rules with the 
potential model results reads: 

/ Bc ~ (2.94 ±0.12)/,, (88) 

which leads to the leptonic decay rate into tv t of about (3.0±0.4) x (V r cb /0.037) 2 x 10 10 s _1 

We have also studied the semileptonic decay of the B c mesons and the g 2 -dependence 
of the form factors. We found that, in all cases, the QCD predictions increase faster 
than the usual pole dominance ones. The g 2 -behaviour of the form factor can be fitted 
with an effective pole mass of about 4.1-4.6 GeV instead of the 6.3 GeV expected from 
a pole model. Basically, we also found that each exclusive channel has almost the same 
rate which is about 1/3 of the leptonic one, a result which is in contradiction with the 
potential model one [Q. Detection of these particles in the next I?-factory machine will 
then serve as a stringent test of the results from the potential models and QSSR analysis. 
The previous analysis is at present extended to the case of the B* meson p7[ . 

10 Conclusion 

We have shortly presented different results from QCD spectral sum rules in the heavy- 
quark sector, which are useful for further theoretical studies and complement the results 
from lattice calculations or/and heavy-quark symmetry. From the experimental point of 
view, QSSR predictions agree with available data, but they also lead to some new features, 
which need to be tested in forthcoming experiments. 



M Bc = 

M B » = 

MA(bcu) = 

MQ( hcs ) = 

M E * (ecu) = 

M s *(bbu) = 
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